Relativity in Fundamental Astoronomy 
Proceedings lAU Symposium No. 261, 2008 
S.A.Klioner, K.P.Seidelman & M.H. Soffel, eds. 



© 2008 International Astronomical Union 
DOI: 00.0000/XOOOOOOOOOOOOOOOX 



Black holes in Active Galactic Nuclei 

M. J.Valtonen\ S. Mikkola^ D. Merritt^, A. Gopakumar 3'®, 
H. J. Lehto^, T. Hyvonen^, H. Rampadarath^'^, R. Saunders^, 
M. Basta^ and R. Hudec^^ 

^Tuorla Observatory, Department of Physics and Astronomy, University of Turku, 21500 

Piikkio, Finland 

^ Centre for Computational Relativity and Gravitation, Rochester Institute of Technology, 78 
Lomb Memorial Drive, Rochester, NY 14623, USA 
■^Tata Institute of Fundamental Research, Mumbai 400005, India 
Joint Institute for VLSI in Europe (JIVE), Postbus 2, 7990 AA Dwingeloo, The Netherlands 
^Leiden Observatory, Leiden University, P.O. Box 9513, NL-2300 RA Leiden, The Netherlands 
^ Department of Physics, University of the West Indies, St. Augustine, Trinidad & Tobago 
Astronomical Institute, Academy of Sciences, Fricova 298, 25165 Ondrejov, Czech Republic 
^Theoretisch-Physikalisches Institut, Friedrich-Schiller-Universitat Jena, 
Max-Wien-Platz 1, 07743 Jena, Germany 
^ Czech Technical University in Prague, Faculty of Electrical Engineering, Technick 2, 166 27 

Praha 6, Czech Republic 

Abstract. Supermassive black holes are common in centers of galaxies. Among the active 
galaxies, quasars are the most extreme, and their black hole masses range as high as to 6- W^^Mq . 
Binary black holes are of special interest but so far OJ287 is the only confirmed case with 
known orbital elements. In OJ287, the binary nature is confirmed by periodic radiation pulses. 
The period is twelve years with two pulses per period. The last four pulses have been correctly 
predicted with the accuracy of few weeks, the latest in 2007 with the accuracy of one day. 
This accuracy is high enough that one may test the higher order terms in the Post Newtonian 
approximation to General Relativity. The precession rate per period is 39°. 1 ± 0°.l, by far the 
largest rate in any known binary, and the (1.83 ± 0.01) • 10^° Mq primary is among the dozen 
biggest black holes known. We will discuss the various Post Newtonian terms and their effect 
on the orbit solution. The over 100 year data base of optical variations in OJ287 puts limits on 
these terms and thus tests the ability of Einstein's General Relativity to describe, for the first 
time, dynamic binary black hole spacetime in the strong field regime. The quadrupole-moment 
contributions to the equations of motion allows us to constrain the 'no-hair' parameter to be 
1.0 ± 0.3 which supports the black hole no-hair theorem within the achievable precision. 
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1. Introduction 

Centers of galaxies typically host dark massive bodies which are thought to be black holes. 
The black hole concept plays a particularly strong role in theories of quasars and other active 
galactic nuclei. Some models use the assumed black hole properties of the dark bodies, such 
as the frame dragging around a spinning black hole (e.g. Camenzind 1989), but on the whole 
the activity is largely explained by accretion onto a central heavy body, as was first pointed 
out by Salpeter (1964) and Zeldovich (1964). Theories of accretion disks have been developed, 
in particular the a disk theory of Shakura and Sunyaev (1973) and its extension to magnetic 
disks by Sakimoto and Corotini (1981) which describe well the origin of activity around massive 
central bodies, whatever the exact nature of these bodies may be. 

In order to prove that the central body is actually a black hole we have to probe the gravita- 
tional field around it. One of the most important characteristics of a black hole is that it must 
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satisfy the so called no- hair theorem or theorems (Israel 1967, 1968, Carter 1970, Hawking 1971, 
1972; see Misner, Thorne and Wheeler 1973). A practical test was suggested by Thorne (1980) 
and Thorne, Price and Macdonald (1986). In this test the quadrupole moment Q of the spinning 
body is measured. If the spin of the body is S and its mass is M, we determine the value of q in 

Q = -9]£- (1-1) 

For black holes q = I, for neutron stars and other possible bosonic structures q > 2. This is 
an important test for stellar mass bodies where stable configurations of various kinds may exist 
(Wex and Kopeikin 1999), but it is also a prime test for the supermassive black hole concept in 
Active Galactic Nuclei and even in our Galactic Center (Will 2008). 



2. OJ287 

BL Lacertae object OJ287 is known to have a quasiperiodic pattern of outbursts at 12 year 
intervals (Sillanpaa et al.(1988),Valtonen et al. (2008c)). Further, the light curve has a double 
peak structure, with the two peaks separated by 1 - 2 years. The available information is that 
the radiation at the peaks is thermal bremsstrahlung radiation, in contrgist to the non-thermal 
synchrotron emission at 'normal' times. The origin of the thermal bremsstrahlung is thought to 
be an impact of a secondary black hole on the accretion disc of the primary (Lehto & Valtonen 

(1996) ). The model has been successful in predicting future outbursts: the predictions for the 
beginning of 1994, 1995 and 2005 outbursts were correct within a few weeks. In addition to 
timing the impacts on the accretion disk, the model includes delays of outbursts relative to the 
disk crossing. The outburst begins when a bubble of gas torn off the accretion disk becomes 
transparent at optical wavelengths. For the 2005 outburst, it also was necessary to consider the 
bending of the accretion disk caused by the secondary (Sundelius et al. (1996),Sundelius et al. 

(1997) ); when combined with the earlier model for the radiation burst delay Lehto & Valtonen 
(1996) the beginning of the 2005 outburst was expected at 2005.74. This is very close to the 
actual starting time of the outburst (Valtonen et al. (2006a), Valtonen et al.(2008a)). It confirmed 
the need for relativistic precession since without the precession the outburst would have been 
a year later. Finally, the gravitational radiation energy loss was included in the prediction for 
the next outburst in September 2007 (Valtonen (2007), Valtonen et al. (2008a)). Observations 
confirmed the correctness of the prediction within the accuracy of one day Valtonen et al. (2008c). 
The model that does not incorporate the gravitational radiation reaction effect is clearly not 
tenable since it predicted the outburst three weeks too late. The probability of so many major 
outbursts happening at the predicted times by chance is negligible. 

Recently, models have been developed which include the spin of the primary black hole (Val- 
tonen et al. 2009). Here we will review this work and discuss its general implications. 



3. Observations 

The model requires accurate timing of the outbursts. The outburst is thought to arise from a 
hot bubble of gas which has cooled to a point where it becomes suddenly optically transparent. 
When such a bubble is viewed from a distance, the emission is seen to grow in a specific way as 
the observational front advances into the bubble. The size of the bubble, and thus also the rate 
of development of the outburst light curve, is a known function of distance from the center of 
the accretion disk. In Figures 1 and 2 we display theoretical light curves for different outbursts 
and overlay them with observations. The labels inside the panels tell the timing which is related 
to the image. 
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Figure 1. The observation of tiie briglitness of OJ287 in four different seasons: 1912/3, 1947, 
1984 and 2005. The heavy curve is a fit through the observed points (dots). Each panel labels 
the deduced start of the outburst. 
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Figure 2. The observation of the brightness of OJ287 in four different seasons: 1972/3, 
1982/3, 1995 and 2007. See caption of Fig. 1. 



4. PN-accurate orbital description 

We invoke the 2.5PN-accuratc orbital dynamics that includes the loading order general rela- 
tivistic, classical spin-orbit and radiation reaction effects for describing the temporal evolution 
of a binary black hole (Kidder (1995)). 
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The 2.5PN-accurate equations of motion can be written schematically as 



X = = *0 + XlPN + XSO + Xq 

+X2PN + X2.r,PN , 



(4.1) 



where x = xi — X2 stands for the center-of-mass relative separation vector between the black 
holes with masses mi and 7712 and xo represents the Newtonian acceleration given by xo = 
— ^-j^ x; m = mi + 7712 and r = \x\. The PN contributions occurring at the conservative IPN, 
2PN and the reactive 2.5PN orders, denoted by iipjv, X2pn and X2.5PN respectively, arc non- 
spin by nature. The explicit expressions for these contributions, suitable for describing the binary 
black hole dynamics, were derived for the first time in the harmonic gauge (Damour (1982)). 
These are well known and they are not repeated here. 

The leading order spin-orbit contributions to x, appearing at 1.5PN order (Barker & O'Connell 
(1975)), reads 



XSO 



(4.2) 



12 [si • (n X v)] 



n + 



(9-|-3v'l-4??) r 



(n X si) — 



7 + Vl-4r7 



{v X Sl) 



where the vectors n and v are defined to be n = x/r and v = dx/dt, respectively, while 
r = dr/dt = n-v, v = \v\ and the symmetric mass ratio ij = mi m2/m^. The Kerr parameter x 
and the unit vector si define the spin of the primary black hole by the relation Si = G m\ x^i/c 
and X is allowed to take values between and 1 in general relativity. Further, the above expression 
for xso implies that the covariant spin supplementary condition is employed to define the 
ccuter-of-mass world line of the spinning compact object in the underlying PN computation 
(Kidder (1995)). Finally, the quadrupole-monopole interaction term xq, entering at the 2PN 
order (Barker & O'Connell (1975)), reads 



2 3 0-^ m?m 

XQ = -qx 2c^/^ 



5(n • Sl 



n — 2(n • si)si 



where the parameter q, whose value is 1 in general relativity, is introduced to test the black hole 
'no-hair' theorem. The precessional motion for the primary black hole spin is dominated by the 
leading order general relativistic spin-orbit coupling and the relevant equation reads 



asi 



n = 



Gmri\ /7-I- VI -4?7 



2c2r2 



(n X v) . 



(4.3) 



It should be noted that the precessional equation for the unit spin vector si enters the binary 
dynamics at IPN order, while the spin contribution enters x at the 1.5PN order. 

The main consequence of including the leading order spin-orbit interactions to the dynamics 
of a binary black hole is that it forces both the binary orbit and the orbital plane to precess. 
Moreover, the orbital angular momentum vector, characterising the orbital plane, precesses 
around the spin of the primary in such a way that the angle between the orbital plane and the 
spin vector si remains almost constant (roughly within ±0.°5 in our model). The spin- vector 
itself precesses drawing a cone with an opening angle of about 12 degrees. 

The precessional period for both the orbital plane and the spin of the binary is about 2400 
years. The orbital inclination relative to the plane of symmetry of the spinning black hole ( and 
presumably relative to the accretion disk of the primary as well) is taken to be 90 degrees. 



5. Timing experiments 

In previous work (e.g. Valtonen 2007) we do the timing experiments using the 1913, 1947, 
1973, 1983, 1984 and 2005 outbursts as fixed points. We cannot determine the Kerr parameter 
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Figure 3. The dots represent solutions witii different values of x (horizontal axis) and outburst 
time (vertical axis). On the y-axis of panel (a) the points scatter around 1995.842, in panel 
(c) around 1982.964 and in panel (d) around 2007.692. In contrast, panel (b) shows a strong 
correlation between the outburst time in 2015 and x- 
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Figure 4. The dots represent solutions with different values of x (horizontal axis) and outburst 
time (vertical axis). On the y-axis of panel (a) the points scatter around 1957.08, in panel (b) 
around 2005.745, in panel (c) around 1912.98 and in panel (d) around 1972.945. 

X and the no-hair parameter q without considering more fixed points. Here we also make use of 
the 1995 and 2007 outbursts which have been observed with high time resolution (see Figure 
and they are therefore suitable for further refinement of the model. In addition, the newly 
discovered 1957 outburst is taken as another fixed point (Valtonen et al. (2006b), Rampadarath 
et al. (2007)). The observed outburst times are listed in Table [1] 



6 



Valtonen et al. 




0.36 



0.36 




0.36 



0.36 



Figure 5. The dots represent solutions witli different values of x (horizontal axis) and outburst 
time (vertical axis). On the y-axis of panel (b) the points scatter around 1959.213, in panel 
(c) around 1984.129, and in panel (d) around 1964.27. Panel (a) displays a strong correlation 
between the timing of the 1945 outburst and X- There are not enough data at present to verify 
the timings of 1959, 1964 and 1945 outbursts, while the 1984 outburst is used as a fixed point 
in our solutions. If new data points are found from years 1945 - 1964, the model can be further 
refined. 



Table 1. Outburst times with estimated uncertainties. These are starting times of the 
outbursts, normalised to the event of 1982.964. 
1912.980 ± 0.020 
1947.283 ± 0.002 
1957.080 ± 0.030 
1972.945 ± 0.012 
1984.130 ± 0.005 
1995.842 ± 0.0015 
2005.745 ± 0.012 
2007.692 ± 0.0015 



With this full set of outbursts no solutions were found unless the parameter x is in the range 
0.2 - 0.36. The majority of solutions cluster around x ~ 0.29, with a one standard deviation 
of about 0.04. The timing experiments give a unique solution for the system parameters: 
precession in the orbit plane per period A(jf>, masses mi and m2 of the two black holes, spin 
parameter x, initial phase 00, initial apocenter eccentricity eo, 'no-hair' parameter q and time- 
delay parameter td- The results are shown in Table [J] 

Here the precession rate is defined as the average change of the apocenter phase angle over 
the last 150 yrs. The time delay parameter depends on the structure of the accretion disk, and 
obtains the value of unity in the model of Lehto & Valtonen (1996). Its value is related to the 
thickness of the accretion disk, as explained in Valtonen et al. (2006b). Figures 3-6 show the 
distribution points, each representing a solution, in four panels per figure. Each panel displays 
the outburst time as a function of X- In Figure 3 the panels (a), (c) and (d) are essentially 
scatter diagrams, demonstrating that the solutions cover the allowed range rather well. On the 
contrary, the panel (b) shows a strong correlation between the time of the 2015 outburst and X- 
The range of possible outburst times extends from early November if x = 0.36 to late January 
2016 if X = 0.2. 
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Figure 6. The dots represent solutions with different values of x (horizontal axis) and outburst 
time (vertical axis). On the y-axis of panel (a) the points scatter around 1971.126, in panel (b) 
around 1934.35 and in panel (d) around 1947.283. In contrast, panel (c) shows a dependence 
of the half-period of the innermost stable orbit on x- The observed half-period is between 43 
and 49 days, implying that x is likely to be at least 0.29. The 1947 outburst is one of our fixed 
points, while there are not enough data to verify 1934 and 1971 outbursts at present. 

Table 2. Solution parameters. 





39.1 ±0.1 


mi 


(1.83 ± 0.01) • 1O^°M0 


m2 


(1.4 ±0.1) • 10* Mo 


X 


0.28 ± 0.08 


00 


56°.5 ±1°.2 


Co 


0.6584 ± 0.001 


1 


1.0 ±0.3 


td 


0.75 ± 0.04 



6. More OJ287's? 

How unique is OJ287? Can wc hope to observe more similar systems? OJ287 is highly variable 
mostly due to strong beaming. At its faintest it goes to ms = 18 which may be taken as its 
intrinsic (unbeamed) brightness. There are about 2 • 10* quasars in the sky brighter than this 
magnitude limit (Arp 1981). Many of these quasars host binary black holes, perhaps as many 
as 50% (Comerford et al. 2009). Thus potentially there are about 10* bright binary quasars 
in the sky to be discovered. But arc they all similar to OJ287? OJ287 is in the last stages of 
inspiral, with only 10* yr left out of its potential 10^ yr lifetime (Volonteri et al. 2009). Thus the 
number of bright short period binaries which could be discovered by techniques similar to the 
discovery of OJ287 is only of the order of ten over the whole sky. Another property of OJ287 
which may be relevant to the discovery rate, is its mass at the upper end of black hole mass 
range. Even though a high mass generally means high luminosity and sampling of large volumes 
of space (observed numbers increasing roughly as M^'^ for luminosity being a constant fraction 
of Eddington luminosity), the frequency function of black hole mass falls with increasing mass. 
The power law slope at the high mass end is about -2.6, i.e. the probability density for mass 
being greater than M has a power law index of about -1.6 (Vestergaard and Osmer 2009). Thus 
the effect of frequency function is to compensate for the volume effect (within the uncertainties). 
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For a binary system with a constant period, there is also another factor. The remaining lifetime 
of the binary goes down as M~'^^^ . The most likely value for a mass from this distribution is 
one quarter of the upper limit (which is of the order of 6 ■ IO'^'^Mq, Vestergaard et al. 2008). It 
is not too far from our measured mass value. The main reason for the discovery of the OJ287 
binary in a sample of about 50 quasars studied by the Tuorla monitoring group is the fact that 
it is a blazar. Because of this, it was a prime target for variability studies. The double peak 
optical outbursts would have easily gone unnoticed if it did not attract attention by its variable 
jet brightness which is seen all the time. Therefore one should extend the periodicity searches to 
quasars which are not highly variable, in hopes of detecting occasional double peaked outbursts. 

From the measured frequency of over 10^" Mq black holes in quasars one easily estimates that 
there must be a dozen or more such quasars in the local volume encompassing OJ287. Since 
their jets would not typically point toward us, they would only appear bright during the short 
intervals of disk impacts ; it would require a major program of historical plate analysis and 
observational monotoring to catch these potentially OJ287-type quasars in action. 



7. Parameter values 

We have already pointed out that the primary mass mi = (1.83 ± 0.01) • lO^^M© is quite 
consistent with what might be expected on general grounds. What about the secondary mass? 
Its value was calculated by Lehto and Valtonen (1996) as m2 « 1.23- 10* Mq, updated to today's 
Hubble constant, while Sundclius et al. (1997) preferred a 30% greater value. These estimates 
are based on the astrophysics of the disk impacts and on the amount of radiation produced in 
these events, and thus they are totally independent of the orbit model. These estimates agree 
well with our more exact value = (1.4 ± 0.1) • 10* Mq. 

The eccentricity of the orbit, if defined by using pericenter and apocenter distances as for a 
Keplerian orbit, is e = 0.7. The initial eccentricity at the beginning of the final inspiral, say, 
10^ yr prior to merger, must have been about e = 0.9, a reasonable value at the initial stage of 
merging black holes (Aarseth 2008). Then by today the eccentricity would have evolved to its 
current value. 

The parameter td gives the ratio of ag to the accretion rate in Eddington units. The value 
td = 0.75 ± 0.04 implies 14 ± 0.5 for this ratio. The mass accretion rate may be about 0.005 of 
the Eddington rate (Bassani et al. 1983) which gives ag = 0.07 ± 0.03, a reasonable value. 

We infer that the primary black hole should spin approximately at one quarter of the maxi- 
mum spin rate allowed in general relativity. There is an additional observation which supports 
this spin value. OJ287 has a basic 46±3 day periodicity (Wu et al.(2006)), which may be related 
to the innermost stable orbit in the accretion disk. However, since we presumably observe the 
accretion disk almost face on, and there is an m = 2 mode wave disturbance in the disk, this 
is likely to refer to one half of the period. Considering also the redshift of the system, and the 
primary mass value given in Table 2 this corresponds to the spin of x = 0.35 ±0.06 (McClintock 
et al. (2006)). The uncertainty of 0.06 units is related to the width of a trough in the structure 
function of the flux variations. For comparison, it has been estimated that x ~ 0.5 for the black 
hole in the Galactic center (Genzel et al. 2003). 

The values of q cluster around q = 1.0 with a standard deviation of 0.1 and maximum range 
±0.3. This is the first indication that the 'no-hair' theorem is valid. Even though there are 
no other known stable configurations than a black hole in the supermassive scale, it is still 
interesting that our result converges at the proper value for general relativity. 
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